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INTRODUCTION and PROBLEM STATE MEMT 

MOTIVATION 


• The OSIRIS-REx mission (Origins Spectral Interpretation Resource 
Identification Security Regolith Explorer) is an asteroid sample 
return mission to Bennu (RQ36) that is scheduled to launch in 
2016. 

• The planned science operations precluding the sample retrieval 
involve operations in terminator orbits (orbit plane is perpendicular 
to the sun). 

• Over longer durations the solar radiation pressure (SRP) perturbs 
the orbit causing it to preccess. 

OUR WORK INVOLVES: 

• Modeling high fidelity SRP model to capture the perturbations during 
attitude changes. 

• Design a stable orbit from the high fidelity models to analyze the stability 
over time. 



INTRODUCTION and PROBLEM STATE MEMT 

PROBLEM STATEMENT 


• Develop High-fidelity SRP model (using SPAD tool) 

• Considers absorption, specular reflection, and diffuse reflection 

• Uses ray-tracing method to model path of sunlight (on OSIRIS-REx 
spacecraft model) 

• Propagate Orbit B in GMAT using: 

• High-fidelity SRP force model (from SPAD) 

• 16 th -order gravity model 

• 3 rd body perturbations: Sun, Jupiter, Earth, 
and Moon 

• Compare with cannonball model 

• Design a stable orbit using Variation 
of parameters (VOP). 




Kyle Hughes, Purdue University 
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SOLAR RADIATION PRESSURE MODELING 



• What is Orbit B: 


OSIRIS-REx's Orbit B 


• Circular 1 km orbit about Bennu (altitude of ~3 R B ) 

• Terminator orbit (orbit normal vector directed at Sun) 

• Key Activity: gather optical data for 12 TAG sites and map gravity field 
• OSIRIS-REx orientation during Orbit B 

• Nominal: nadir pointing, with solar panels directed at Sun 

• HGA Pointing: Each (Earth) day, HGA pointed at Earth (for 8 hrs), with body-fixed 
z axis pointed at Sun (approximately) 

• Attitude changes cause change in SRP force on spacecraft — which needs to be 
modeled for a high-fidelity trajectory simulation 
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SOLAR RADIATION PRESSURE MODELING 

Solar Radiation Force Model 


?SRP ~?a+ ?s + ?d 

= -P srp Acos0 {(1 - C 5 )S + 2 (Cj-cosS + ;Cd)w] 
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• Ray tracing is the process of following 
the path of a light ray as it absorbs 
and/or reflects off of various surfaces of 
an object. 

• Used to model the SRP forces acting on 
the spacecraft in high fidelity 

• Create pixel array normal to Sun- 
spacecraft vector 

• Project points at each pixel center onto 
spacecraft surface 

• Compute SRP force at each surface 
contact point and find net force 



Ziebart, M., and Dare, P., “Analytical Solar Radiation Pressure Modeling 
for GLONASS using a Pixel Array”, Journal of Geodesy, Vol. 75, 2001 . 


SOLAR RADIATION PRESSURE MODELING 



OSIRIS-REx CAD Model for SPAD Analysis 


• Ray tracing modeled using the Solar 
Pressure and Aerodynamic Drag 
(SPAD) tool, developed at GSFC in 
Code 595, Navigation and Mission 
Design Branch 

• The tool takes in inputs of: 

• Material properties (emissivity etc) 

• Dimensions (to generate CAD model) 

• Pixel size (for light-ray grid) 

• Azimuth and elevation (S/C model SPAD-generated OSIRIS-REx model 

orientation) 

• SPAD outputs the Force/SRP 

• Trajectory propagation tool (GMAT) 
computes force on S/C using computed 
Sun-S/C distance, and specified S/C 
attitude. 
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VARIATION OF PARAMETERS 

Variation of parameters 


Assumptions 

• The following assumptions will be made for the purposes of this 
analysis. 

• The orbit of the spacecraft about the asteroid is high enough that an 
inverse square law asteroid gravity model suffices. (For the case of 
OSIFtlS-Rex Orbit B, the orbital radius is roughly equal to four 
asteroid radii, so this assumption is certainly valid.) 

• The eccentricity of the spacecraft orbit about the asteroid is small 
enough that first-order modeling of eccentricity effects is valid. 

• The spacecraft is in a near-terminator orbit about the asteroid. 

• The spacecraft has a constant ratio of projected area to the Sun and 
mass, and behaves as a pure solar absorber. 

• The changes in the spacecraft orbital elements that are caused by 
solar radiation pressure occur slowly enough that integration over a 
complete spacecraft orbit, using constant element values in the VOP 
equations across this range, is valid. 



VOP Equations 


VARIATION OF PARAMETERS 
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VARIATION OF PARAMETERS 

AVERAGED PER-REV ELEMENT CHANGES 


• The components of SRP perturbation acceleration in the 
spacecraft rotating RSW frame can be computed and put into 
the Gaussian VOP equations, giving the instantaneous 
variations in each orbital element 

• Since the spacecraft orbit has e « 1 , linearization of 
eccentricity effects is acceptable 

• Integrating over a complete spacecraft rev gives the per-rev 
changes in the spacecraft orbital elements: 

a ~ 0, 

3 ( / )(1 e 2 ) 3 a 2 p s (c r A / m) sin sin , 

i~3( / )e(l e 2 ) 2 a 2 p s (c r A / m) cos cos , 

«3( / )e(l e 2 fa 2 p s (c r A / m) cosec z'cos sin , 

«3( / )e 1 (1 e 2 ) 2 a 2 p s (c r A/m ){ sin cos e 2 cot i cos sin }. 



VARIATION OF PARAMETERS 

USE FOR STABLE TERMINATOR ORBIT DESIGN 


• Desired geometry of spacecraft orbit about asteroid: 

• Angle A either 0 (giving CCW spacecraft motion about asteroid as 
viewed from the Sun) or 180 deg (CW motion) 

• Inclination (relative to orbital plane of asteroid) 90 deg 

• All elements constant except for RAAN: this must precess at 
asteroid orbital rate, in order for orbit to remain on terminator 

• The conditions on A and i imply that the per-rev changes in e 
and co will be zero, as desired 

• Zero per-rev change in / is obtained by setting co to either 90 
deg (CCW motion case) or 270 deg (CW) 

• Selecting the correct eccentricity value (see next slide for 
expression) can then be used to adjust the per-rev change in 
/2to give the desired nodal regression rate 

• Mechanism: if orbit is eccentric (here, mildly so), SRP gives a 
net out-of-plane Av centered at apogee: this torques orbit, 
changing /2(but not i, for ryas selected) 
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RESULTING ORBIT DESIGNS 

ORBIT B SCENARIO 


SRP-Stable Orbit B 

• A SRP-stable orbital eccentricity value of 0.139, using the following 
numerical values that correspond to OSIRIS-Rex Orbit B: 

• Spacecraft 

• mass (during Orbit B): m = 1 ,198 kg 

• projected area: A = 1 2 m 2 

• coefficient of reflectivity: c r = 1.4 

• mean orbital speed: v orb = 0.06450 m/s 

• Bennu gravitational parameter: p = 4.16 km 3 /s 2 

• Spacecraft-orbit semi-major axis (about Bennu): a = 1 ,000 m (4 times 
asteroid radius of -250 m) 

• Bennu’s mean motion (heliocentric): n= 1.67x10' 7 rad/s 

• Bennu’s solar distance: 1.11 AU 

• SRP:p s = 3.69x1 O' 6 N/m 2 

• The corresponding spacecraft apoapse and periapse radii are 1,139 m, 
and 861 m, respectively. 
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Cannonball 


RESULTING ORBIT DESIGNS 



• In this view, an un-perturbed terminator 
orbit would appear as a single circle or 
ellipse. The yellow arcs shown in the 
figures represent times when the 
spacecraft is in its nadir-pointing 
orientation (16 hrs per day), and the 
blue arcs represent times when the 
spacecraft is oriented with its HGA 
directed at Earth (8 hrs per day). For a 
spherical spacecraft however, the 
orientation changes make no impact on 
the SRP force or the trajectory. 

• Figure A shows the propagated Orbit B 
trajectory using the nominal initial orbit 
conditions (with e = 0), and Figure B 
shows the Orbit B propagation with the 
redesigned orbit (e = = 0.139). Over the 
60-day propagation, both trajectories 
exhibit perturbation; however, the 
redesigned trajectory in Figure B is 
clearly improved. 



Orbit B Propagations: Cannonball SRP Model, A = 12 m 2 
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RESULTING ORBIT DESIGNS 


SPAD MODEL 


The top figures show the trajectory in the 
Bennu-Sun rotating frame (similar to that of 
Figure 6), where the view is along the Bennu- 
Sun direction. The propagation of the Orbit B 
trajectory using the nominal initial orbit 
conditions (with e = 0) is shown in Figure A and 
the propagation of the stabilized Orbit B (with e 
= = 0.139) is shown in Figure B. 



Both trajectories are clearly perturbed; however, 
the designed trajectory in Figure B exhibits less 
perturbation over the 60-day simulation. The 
yellow arcs shown in the figures represent 
times when the spacecraft is in its nadir- 
pointing orientation, and the blue arcs represent 
times when the spacecraft is oriented with its 
HGA directed at Earth. Unlike with the 
cannonball SRP model, however, the 
orientation changes do have an impact on the 
SRP force, and thus, the resulting propagation 
of Orbit B. 


Orbit B Propagations: SPAD SRP Model 





Conclusions 


/ JT 

QSIRIS-R^K \ 

• SRP plays a significant role in maintaining the nominal 
terminator orbit about Bennu 

• Significant differences in trajectories found between 
cannonball and SPAD SRP models 

• Designed new SRP-stable terminator orbit using VOP 
equations, which requires an eccentricity of about 0.139 
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The results show that there are significant differences between the three 
trajectories investigated for this analysis — one key difference being the variation 
in eccentricity. 

The SPAD and 12-m2 cannonball cases exhibit a lag in Q precession (compared 
to the 20-m2 cannonball case). A key contributor to this precession lag seems to 
be the smaller amount of SRP force experienced by the SPAD and 12-m2 
cannonball models (which is due to their smaller cross-sectional areas). 

The lag in Q precession increases the deviation between the orbit normal and 
the Sun direction, which is shown by the increase in the amplitude of |3 (the 
angle between the orbit normal vector and the Sun line). 

Larger values of |3 cause more of the SRP force to be directed in the orbit- 
relative radial and transverse directions. The transverse force, 
the change in eccentricity. 

FSRP, T 
FSRP , N 

FSRP, T raises apoapsis (increasing e, increasing a) 

FSRP, T lowers periapsis (increasing e, decreasing a) 




Influence of SRP on OSIRIS-REx at Bennu 


Acceleration due to gravity of Bennu and Sun: (in nominal Orbit B) 
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Accel due to Bennu’s gravity * 3.88 x 10 6 m/s 2 
Accel due to Sun’s gravity * 3.93 x 10' 11 m/s 2 
Net Accel due to gravity * 3.88 x 1 0 6 m/s 2 


• Acceleration due to SRP (~1% the accel due to gravity) 

Solar pressure at 1 .01 5 AU, P SRP = 4.49x1 O' 6 m/s 2 mass of SC = 1198 kg 
Cannonball (A=12m 2 , C R =1.4): a can =C R A P SRP / m ~ 6.30 x 10‘ 8 m/s 2 

SPAD Model (Nadir): a Nadir * 5.29 x 10' 8 m/s 2 

SPAD Model (Earth-Pointing): ^Earth-Pointing * 5.1 x 10‘ 8 m/s 2 
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